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We present an efficient ab initio approach for the study of magnetic transport properties based 
on the Boltzmann equation with the Wannier interpolation scheme. Within the relaxation time 
approximation, band-resolved electric conductivity under a finite magnetic field is obtained and 
the historical motion of the electron wave packet in reciprocal space is determined. As a typical 
application of this method, we have calculated the electric conductivities of MgB2 under finite 
magnetic fields. Multiband characters for the individual bands are revealed, and the field dependence 
of the conductivity tensor is studied systematically with the field orientated parallel and normal 
to the c-axis, respectively. The obtained historical motion is employed to simulate directly the 
cyclotron motion in the extremal orbit and determine the corresponding effective mass. Moreover, 
This approach is further exploited to calculate the Hall coefficient in the low-field limit, without the 
complicated computation for the second k derivative of the band. 



PACS numbers: 71.15.Dx, 71.18.+y, 75.47.-m 

I. INTRODUCTION 

Magnetic transport properties, such as magnetoresis- 
tance (magnetoconductivity) and the Hall effect, which 
have been investigated experimentally and theoretically 
for a long history, are of great importance in the field 
of research for metals and semiconductors. These prop- 
erties are sensitively dependent on the topology of the 
Fermi surface (FS). The corresponding semiclassical in- 
terpretations were first given by Lifshitz et aZP in the 
light of the Boltzmann equation with a general collision 
operator. Subsequently they we re f urther developed for 
the studies of magnetoresistance^' and Hall effect,^ re- 
spectively. 

The Boltzmann transport theory has functioned suc- 
cessfully in previous studies for two-dimensional systems 
through analytical and numerical solutions Yet in real 
three-dimensional crystals, transport behaviors are sig- 
nificantly determined by the complex FS, where models 
are not much satisfactory and analytical solutions are 
formidable. Ab initio study is then essential to provide 
the most comprehensive and reliable investigation. To 
our knowledge of magnetoresistance, so far works based 
on electronic structure calculation have been rarely con- 
cerned with the numerical solution of the Boltzmann 
equation for the band-resolved electric conductivity un- 
der the finite magnetic field I21i0j They were faced with 
either heavy ab initio computation of properties on the 
accurate FS, or unreliable results from the Shankland- 
Koelling-Wood (SKW) interpolation scheme^^ j n the 
presence of band crossings close to the FS. On the other 
hand, in the low-field limit the solution of the Boltzmann 
equation has been extensively studied for the Hall ef- 
fect through the power-series expansion of the magnetic 
fieldP^ With the isotropic relaxation time approxima- 
tion, the Hall coefficient is explicitly formulated for met- 
als having cubic symmetry. 4 Following this perturbation 
method, previous simulations were carried out based on 
empirical tight-binding (TB) band structures,^ ab initio 



calculation combined with the SKWp^ and the Wannier 
interpolation scheme!^ However, only cubic metals were 
involved in these studies, and the second k derivative of 
the band was necessary. 

In this paper, we present an ab initio approach com- 
bined with the Wannier interpolation scheme for the 
solution of the Boltzmann equation within the relax- 
ation time approximation. Band structures are obtained 
through the ab initio calculation, and Wannier functions 
(WFs) are then constructed.^ The Runge-Kutta scheme 
is employed to determine the historical motion of the 
electron wave packet, where the WFs are used to evalu- 
ate the local velocities. The numerical integration with 
the tetrahedron scheme^ is adopted to deal with the sum 
over k-space. It is noticeable that the velocity evaluation 
with the aid of the Wannier interpolation is more efficient 
(reliable) than that using ab initio calculation ( the SKW 
interpolatiorf^J) at individual k points. We have already 
applied this magnetoconductivity simulation to the su- 
perconductor MgB 2 ^ in our previous Letter^ and the 
combination with the experimental measurements played 
a key role in deriving electron-phonon scattering times. 
In this work, we further study systematically the band- 
resolved electric conductivity of MgB2 based on this ap- 
proach, with the finite magnetic field orientated parallel 
and normal to the c-axis, respectively. The historical 
tracks obtained during the Runge-Kutta integration is 
also used to study directly the cyclotron motion of the 
electron wave packet in the extremal orbit. This provides 
a direct approach to determine the effective mass of the 
corresponding extremal orbit. Moreover, we are able to 
calculate the Hall coefficient in the low-field limit based 
on the conductivity simulation, which is applicable for 
crystals of arbitrary structures. Unlike the perturbation 
method, the complicated computation of the second k 
derivative of the band is avoided. 

The paper is organized as follows. The simulation 
method is presented in Sec. II. It is employed to study 
systematically the magnetoconductivity of MgB2 in the 
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next section, where the simulation of the cyclotron mo- 
tion of the electron wave packet is included. In Sec. IV, 
we calculate the low-field Hall coefficients for several cu- 
bic metals and hep Mg. Conclusions are drawn in section 
V. 



II. SIMULATION METHOD 



With the powerful Wannier interpolation scheme ES122I 
we calculate the band-resolved electric conductivity un- 
der a finite magnetic field through numerically solving 
the Boltzmann equation. In this section, all the ingredi- 
ents involved in our simulation are detailed below. 



A. Band-resolved electric conductivity 

According to the Boltzmann transport theory within 
the relaxation time approximation, the dc electric con- 
ductivity tensor of a metal under a uniform magnetic 
field is written as^ 
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where n is the band index, r n (e n (k)) is the scattering 
time of the n-th band and assumed to be dependent on 
the wave vector k only through the eigenvalue e„(k). / is 
the Fermi-Dirac distribution. v„(k) is the velocity given 
by the equation of motion 
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v n (k) is a weighted average of the velocity over the past 
history of the electron passing through k: 
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The magnetic field is explicitly acting on the time evolu- 
tion of k n (f), as given in the equation of motion: 



dk n (f) e 
df = -jjV n (K(t)) x B . 



(4) 



with the initial condition k n (0) = k. We start with this 
condition for one k point and simulate reversely the time 
evolution given by Eq. (4). It is physically reasonable 
and mathematically maneuverable, and follows a typical 
Runge-Kutta method to deal with the equations of mo- 
tion. The velocity at an arbitrary k points is evaluated 
analytically through the Wannier interpolation™ as de- 
tailed in the following subsections. A serial of k(t) points 
is obtained as a result of the Runge-Kutta integration for 
each initial k point. They sample the history track of the 
electron wave packet at the corresponding time. Conse- 
quently, v(k), according to Eq. (3), can be obtained as 



the weighted average of the velocities at these k(f) points 
along the track. Then the integration in Eq. (1) is per- 
formed by employing the modified tetrahedron scheme.^ 

As in most cases, the relaxation time for each band 
is not easily determined. Therefore, we incorporate the 
magnetic field B and the relaxation time r as one variable 
Bt, which relates to an important dimensionless quantity 
lot = ^prBr, with m* the effective mass. Phenomenally, 
lot represents the completed part in the cyclotron orbit 
before the carrier is scattered out. 

An implicit rule holds in the semiclassical framework: 
the magnetic field does not affect the kinetic energy of 
the carrier. Thus it imposes a severe restriction on ev- 
ery point in the k(f) serials that they are all sampled 
from the same constant-energy surface. It is a straight- 
forward consequence from the theoretical viewpoint, but 
not accomplished automatically in the numerical simula- 
tion. In our code, we developed a modified self-adaptive 
Runge-Kutta scheme, which takes much less time than 
the regular Runge-Kutta methods to achieve the same 
accuracy. We did convincing tests for several arbitrary 
initial k points and found that the energies in their k(f) 
serials were converged better than Hartree. 



B. Wannier interpolation 

In our simulation, the velocity at an arbitrary k 
point is evaluated through analytical calculation with 
the Wannier interpolation scheme. WFs are constructed 
as a "postprocessing" operation based on Bloch eigen- 
states and eigenvalues obtained through a standard ab 
initio electronic structure calculation carried out on a 
uniform fc-point grid!^ Here we employ the ultrasoft- 
pseudopotential^ plane-wave method with generalized- 
gradient approximation^ for the exchange and correla- 
tion potential. We give a brief review about the Wannier 
interpolation as follows, while details can be found in 
Ref. [17l We use the symbol q to denote the points on 
the ab initio mesh, and k for arbitrary and interpolation- 
grid points. 

The eigenstate \tpnq) is obtained after the electronic 
structure calculation. Its periodic part is defined as 
ft nq (r) = e~ lq ' r ?/! nq (r). Denoting the Wannier function 
in the cell R within the band n as |Rn), we have 
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(5b) 
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Here V is the volumn of the real-space primitive cell. 
As shown in Eq. (5a), the freedom in the choice of the 
phases of the Bloch states results in the arbitrariness of 
the Wannier functions. By minimizing their total delocal- 
ization with respect to the phase freedom, the unique set 
of maximally localized Wannier functions is determined. 
The corresponding Bloch-like functions ) are then 
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obtained according to Eq. (5b). The superscript (W) or 
(H) marks the quantities belon ging t o the Wannier or 
Hamiltonian gauge, respectively! 16 * 20 ! Straightforwardly 
we obtain the Hamiltonian by 

ijW(q) = «)|ff(q)|uW). (6) 

With a Fourier transform, we have 

(On\H\Rm) = -L £ e -^ R ffW(q) , (7) 

« q 

where -ff is the effective one-particle Hamiltonian. The 
sum runs over all the q points on the ab initio mesh. 
Finally the interpolation is implemented by transforming 
back the matrix (On\H\Rm) to an arbitrary k point, 

/^>(k) = £e* k - R <On|iT|Rm). ( g ) 

R 

The Hamiltonian at any k point is then analytically de- 
termined. 




FIG. 1: (color online) The calculated Fermi surface of MgB2 
and the extremal orbits labelled for the magnetic field along 
the c-axis. 



TABLE I: The asymptotic behaviors of the diagonal elements 
of the electric conductivity tensor of MgB2 in the classical 
high-field limit. The coordinate system is defined with the 
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direction. 
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B || x 
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C. Evaluation of velocity 

To calculate the velocity, we diagonalize the above in- 
terpolated Hamiltonian by finding a unitary rotation ma- 
trix U(k): 



H (H \k) = [ft(k) flW(k) t/(k) , 



(9) 



where ff4m (k) 



e (H) S 



and e„ k are identical to the 



eigenvalues obtained from the ab initio calculation at k = 
q. We then evaluate the velocities for an arbitrary k 
point through 
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Note that the WFs are constructed via discrete Fourier 
transform based on the ab initio band structure. The 
interpolation will not be je opardized by band crossings 
and avoided crossings! 16 1 20 1 



III. MAGNETO CONDUCTIVITY UNDER 
FINITE MAGNETIC FIELDS 

As the first remarkable application of our approach, we 
study the electric conductivities of normal-state MgB 2 
under finite magnetic fields. This superconductor has 
been investigated extensively in recent yearsP^HUD i n our 
previous Letter^, the simulation method presented in Sec. 
II was employed and the combination with the experi- 
mental measurements provided a unique way to obtain 
band- resolved scattering rates in multiband systems. In 
this work, we further study thoroughly the anisotropic 
magnetoconductivity of MgB 2 with the magnetic field 
orientated parallel and normal to the c-axis, respectively. 



The FS of MgB 2 consists of four sheets. Two forming 
slightly warped cylinders come from a couple of "hole- 
like" quasi-two-dimensional bands (bonding o\ and an- 
tibonding 02), and two tubular sheets form the three- 
dimensional 7T bands (antibonding 7Ti and bonding 7^), 
as shown in Fig. 1. With the magnetic field along the c- 
axis, there are six extremal orbits centered around T and 
A, as labelled in Fig. 1 by numbers 1-3 and 4-6, respec- 
tively. The coordinate system is set with the c-axis along 
the z direction throughout the conductivity calculation. 

Before calculating the magnetoconductivity of MgB2, 
we first analyze its asymptotic behavior for each band in 
the classical high-field limit. According to the Boltzmann 
transport theory, the change in electric conductivity un- 
der a finite magnetic field arises from the cyclotron mo- 
tion of the electron driven by the field. In general, contri- 
butions from the cyclotron motions in a closed orbit and 
an open orbit have different field dependences. Follow- 
ing the analytical deduction by Lifshitz et aiP, and the 
qualitative analysis by PipparcP^l and FawcettpS based 
on the geometric features of the FS, we can determine 
the asymptotic behavior of the diagonal element of the 
electric conductivity tensor of MgB 2 , as summarized in 
Table 1. It is helpful to check the validity of the following 
numerical simulation under finite magnetic fields. 



A. Magnetic field along the c-axis 

When the field is applied in the c direction, all the 
four Fermi sheets are composed of closed orbits. We use 
Eq. (1) to calculate the electric conductivity along the x- 
axis in the ab plane. o~ xx /t of the four bands are shown 
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FIG. 2: (color online) a xx /r as a function of Bt of the four 
bands of MgB2, with the magnetic field along the c-axis. 



in Fig. 2. 

The simulation is carried out with the FS refined after 
the conventional electronic structure calculation.^ With 
the magnetic field along the 6-fold rotation axis, the sym- 
metry of the system are reduced by half, and a yy , which 
are not plotted here, behave the same as a xx The field 
dependence of the band-resolved a xx can be described 
well with a dominant term B~ 2 ®m good agreement with 
the asymptotic behavior in Table 1. Since the magni- 
tudes of a xx /r of the two ir bands in Fig. 2 are noticeably 
different, understanding of generalized a and n bands 
might be incomplete, and the individual bands deserve 
further studies. 

As stressed in Sec. II, the historical motion of the elec- 
tron wave packet in reciprocal space is obtained in the 
process of the conductivity calculation. We take advan- 
tage of these cyclotron motions in the extremal orbits of 
the four bands to determine the corresponding effective 
masses. It provides a straightforward approach to in- 
vestigate the extremal orbit and effective mass, which is 
usually a major concern of the study of the de Hass-van 
Alphen (dHvA) effect Pt22l 

In one of these previous studies of the dHvA effect 
in MgB2, the calculation was carried out with m* = 
T^r § | v ^| , and provided the absolute values of the effec- 
tive massesPS Later, Mazin et aZP^ obtained the effective 
mass through the standard formula: m* = §^(§g)k, , 
with A the area of the orbit, which was computed in 
small energy ranges around the Fermi level. 

In our case, the cyclotron motion of the electron wave 
packet at the FS is tracked accurately with the serial 
of k(i) points obtained through the self-adaptive Rungc- 
Kutta integration for every initial k point. Particularly, 
in the extremal orbits, as labelled in Fig. 1, we choose 
one initial k point each to study the historical motion. 
As an example, we plot in Fig. 3 two of these historical 
tracks in Orbit 1 and Orbit 3, respectively. 

These sinusoidal-like curves shown in Fig. 3 illustrate 
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FIG. 3: (color online) Historical track of the electron 
wave packet in reciprocal space. The initial positions are 
k(0.065, 0.0, 0.0) in Orbit 1 (a), and k(0.444, 0.0, 0.0) in Orbit 
3 (b), respectively, k points are expressed as k(fci, &2, £3), in 
the unit of reciprocal lattice vectors. There are 4000 points on 
each curve. £3 is fixed in the Y plane during the simulation. 
The abscissa ui £ t is related to Bt as ui e t = ^—Bt, where m e 
is the mass of a free electron. 



TABLE II: Calculated effective masses in the extremal orbits 
of MgB2, given in the unit of the free-electron mass. 



Orbit: 


1 


2 


3 


4 


5 


6 


Ref. ET] 


0.254 


0.313 


1.699 


0.550 


0.612 


0.924 


Ref. El 


-0.251 


-0.543 


1.96 


-0.312 


-0.618 


-1.00 


This work 


-0.255 


-0.541 


1.816 


-0.309 


-0.656 


-0.988 



vividly the cyclotron motion in reciprocal space driven by 
the magnetic field. The period of the motion, T, can be 
read directly from the figure. The effective mass of the 
corresponding extremal orbit is then straightforwardly 
obtained in the absolute value as \m*\/m e = T/2ir, 
where m e is the free-electron mass. Then the sign of 
the effective mass is determined through analyzing the 
direction of the cyclotron motion. As shown in Fig. 3, 
the cyclotron motion in Orbit 3 is anticlockwise, while it 
is clockwise in Orbit 1, as well as those in the rest four 
extremal orbits. It is then concluded that the third or- 
bit is "electronlike", while the rest of all are "holelike". 
Our results, listed in Table 2, show good agreement with 
the previous studies, which confirms the validity of the 
approach. 

In the rest of this subsection, we return to the inves- 
tigation of band-resolved electric conductivity. In the 
current geometry, the Hall conductivity, a xy , is another 
non-trivial element besides a xx . In our previous Letter jSI 
the Hall conductivity was calculated with Eq. (1) for all 
the four bands. A sign change from positive to negative 
has been found in the 7r 2 band with increasing field, while 
the 7ti and the two a bands behave simply as holes. For a 
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FIG. 4: (color online) The partial Hall conductivity 
(Txy (fe z )/r'™- ) contributed by orbits between — fc z to k z of the 
individual bands as a function of k z . 



thorough understanding, we further perform the partial 
integration for the Hall conductivity, a x n y {k z )/T^ n \ by 



4^)/r (n) 



4?r 3 / dk x dk y 



d^")(k)4«)(k) 



(ii) 



We employ a 30 x 30 x 40 mesh, which gives 3192 k 
points in the irreducible Brillouin zone. The extra finer 
grid along the z-axis is adopted for the integration from 
-k z to k z . Bt is fixed at 1.0 x KT 12 Ts. 

The calculated partial Hall conductivity of the four 
bands is shown in Fig. 4. For the two a bands, the lin- 
ear growth is understood with the cylindrical shape of 
the Fermi sheets. For the two ir bands, there are sign 
changes along the k z axis, which confirm that multiband 
character exists in the individual tt bands. Although 
the 7Ti band was considered "electronlike" according to 
the band structure, "holelike" part occupies most of the 
Fermi sheet except for the small peaks close to the T 
plane. Thus it behaves as a hole in the full Hall conduc- 
tivity with the increasing magnetic fieldP Yet while the 
situation is opposite for the tt2 band. The small hills are 
"holclike" , and the big belly is "electronlike" . Thanks to 
the ab initio simulation of band-resolved conductivity, we 
have obtained such detailed knowledge of the multiband 
character for the individual bands in MgB2. 



B. Magnetic field normal to c-axis 

When the magnetic field is normal to c-axis, different 
transport behaviors are expected due to the hexagonal 
structure of MgB 2 . We further study this geometry by 
applying the field along the x-axis in the ab plane. The 
one-dimensional infinity of the two a bands is then ex- 
posed, and the tubular 7r bands consist of closed orbits in 
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FIG. 5: (color online) Bt dependences of g xx /t (a), (Jyy/r 
(b), and a zz /r (c) of the four bands of MgB2. The magnetic 
field is applied along the x-axis in the ab plane. 



the planes normal to the field, including extended closed 
orbits. 

We calculate the electric conductivity tensor using 
Eq. (1). Symmetry is further reduced because the 6- 
fold rotation axis has vanished in the current geometry 
The three diagonal elements of the conductivity tensor 
are shown in Fig. 5. 

There are two categories of field-dependent behaviors 
of the diagonal elements. For the two a bands, conduc- 
tivity elements show little dependence on the magnetic 
field. In the x direction, the motion of the electron wave 
packet is not affected by the field, rendering a xx field- 
independent. Similar behavior in the y direction is ac- 
counted for b y open orbits in the two (nearly standard) 
cylinders 111241 The z-components of Fermi velocities are 
negligible compared to those in the other two directions. 
This provides an elegant interpretation for the almost 
vanished elements o~ zz . The two n bands are composed 
of closed orbits with the magnetic field along the x-axis. 
As the field increases, o~ xx of the two 7r bands decrease 
and tend to saturation in the high-filed limit, while a yy , 
as well as a zz , shows a decay which can be described well 
by a dominant term B~ 2 . All the results agree with the 
high-field asymptotic behaviors listed Table 1. In addi- 
tion, these conductivity elements under finite fields can 
be fitted to obtain global field dependences, as carried 
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out in Ref. M 

In principle, there are two applications of the calcu- 
lated band-resolved electric conductivity under the finite 
magnetic fields. With the determined relaxation time in 
each band, magnetoresistance study can be performed.^ 
Alternatively, combining the calculation results with the 
experimental measurements of conductance or resistance 
in different temperatures can lead to quantitative estima- 
tion about these scattering rates, including the impurity 
scattering and the electron-phonon couplingP It is worth 
pointing out that the efficient calculation provides reli- 
able results and the simulation is performed at the accu- 
rate Fermi level. We suggest that this ah initio approach 
is helpful in the study of finite-field magnetic transport 
properties. 



IV. LOW-FIELD HALL COEFFICIENT 
CALCULATION 

Beyond the electric conductivity simulation under a fi- 
nite magnetic field, we further applied the method in the 
low-field limit, where Hall effect is the major interest of 
the magnetic transport research. Previous studiefSKI] 
followed the perturbation method, which adopts the 
power-series expansion of the magnetic field and gives 
the Hall coefficient for cubic crystals as 
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FIG. 6: The band structure of fee Pd. The Fermi level is fixed 
at eV. 



Here <7o is the conductivity: 
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and an is the Hall conductivity: 
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M" 1 is the inverse mass tensor, and its elements are de- 
fined by 



m; 



1 8 2 e 
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(15) 



In this perturbation method, the evaluation of the sec- 
ond k derivative of the band is a difficult task. It either 
gives unreliable results with the SKW interpolation in the 
presence of band crossings and/or near degeneracies, 15 
or follows a very complicated formula with the Wannier 
interpolation!^ For instance, Fig. 6 shows the complex 
band structure of fee Pd. It has three energy bands 
across the Fermi energy, and multiple band crossings in 
the proximity of the Fermi level. Expensive computation 
is inevitable during the SKW interpolation due to the 



complexity of the band structure, and the accuracy of 
the result will be threatened because of the band cross- 
ings. Even with the Wannier interpolation, the heavy 
computation of the inverse mass tensor still can not be 
avoided. 

In our case, each element of the electric conductivity 
tensor can be obtained through Eq. (1). We can further 
evaluate the Hall resistivity according to the definition 
directlyP2 
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' xx^yy 
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Here we assume that the magnetic field is applied along 
the z-axis. The low-field limit is achieved with p yx grow- 
ing linearly with the field. The Hall coefficient is then 
given by the definition, as wellP 
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(17) 



This scheme is valid for metals with arbitrary structures. 
Particularly, in a crystal with threefold or fourfold sym- 
metry about the magnetic field direction, rotation trans- 
formations require a xx = a yy and a xy — —cr yx . Thus we 
have 



Pyx 



(18) 
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TABLE III: The calculated and experimentally observed Hall 
coefficient Rh, in the unit of 10~ 11 m 3 /C. Experimental data 
are all from Ref . HI 





Li 


Al 


Cu 


Pd 


Mg 


Ref.[H 


-12.8 


-1.7 


-5.2 


-6.0 




Ref. US] 


-12.4 


-3.4 


-4.9 


-17 




Ref.ua 


-12.7 


-2.5 


-4.9 


-11.9 




This work 


-13.0 


-2.2 


-5.2 


-11.2 


-10.4 


Experiment 


-15.5 


-3.43 


-5.17 


-7.6 


-8.3 



to simplify the calculation in the application to all cubic 
and hexagonal systems. Throughout the calculation for 
the low-held Hall coefficient, we do not need to deal with 
the inverse mass tensor. The complexity of the band 
structure will not result in expensive computation. 

We performed the calculation of the low-held Hall co- 
efficient with Eqns. (16), (17), and (18) for typical cubic 
metals Li, Al, Cu, and Pd, as well as hep Mg, where 
the magnetic field is along the hexagonal axis. The re- 
sults are listed in Table 3. For Li, Al, and Cu, we find 



great agreement with previous simulations based on band 
structures fitted through TB schemeP ab initio calcula- 
tion combined with SKW interpolation^ and Wannier 
interpolation.^ With a certain ratio between the relax- 
ation times of the two bands of aP, we can reproduce 
the experimental value in our simulation. In the case of 
Pd, the result with the SKW interpolation is not as sat- 
isfactory as those through TB and with the Wannier in- 
terpolation. It results from the complex band structure 
as analyzed above. Our result agrees excellently with 
that from Ref. [THl yet without the complicated evalua- 
tion of the inverse mass tensor. The discrepancy between 
calculated and experimental values was acc ounte d for by 
the constant relaxation time approximation!^^ For hep 
Mg, our result is comparable to the experimental data 
with the magnetic field parallel to the hexagonal axis. 



Wannier interpolation scheme is employed to analytically 
evaluate the velocities along the historical track of the 
electron wave packet. The simulation is guaranteed to 
be performed at the accurate FS by the modified self- 
adaptive Runge-Kutta method. As the first application, 
we have studied systematically the electric conductivity 
of normal-state MgB 2 with the magnetic field orientated 
parallel and normal to the c-axis, respectively. The multi- 
band characters within the individual bands are then re- 
vealed in detail. 

During the conductivity calculation, the cyclotron mo- 
tion of the wave packet driven by the magnetic field is 
simulated. We have taken advantage of these motions in 
the extremal orbits of MgB 2 to calculate the correspond- 
ing effective masses, with the results in good agreement 
with previous studies. It provides a direct simulation 
of the cyclotron motion and a reliable evaluation of the 
effective mass of the extremal orbit. 

Another application of this method is the calculation 
of the low-field Hall coefficient as a universal approach for 
metals and semiconductors of arbitrary structures. We 
have reproduced the Hall coefficients for several typical 
cubic metals, and for hep Mg as an example of non-cubic 
systems. 

A remarkable appeal of this approach is that it pro- 
vides an ab initio method to simulate the electric con- 
ductivity under a finite magnetic field. It is efficient 
and reliable owning to the Wannier interpolation scheme. 
Moreover, the direct simulation of the cyclotron motion 
is an elegant and useful application. The calculation of 
the low-field Hall coefficient is universal for crystals of 
arbitrary structures, and avoids all complicated compu- 
tation of the second derivatives of the band. Therefore, 
this approach can be generally applied to a wide variety 
of systems for magnetic transport study. 
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